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INTRODUCTION 


This  work  is  part  of  an  investigation  with  the  long-range  objective 
of  predicting  projectile  trajectories  and  the  forces  applied  to  projec- 
tiles during  impacts  against  thin  plates  of  moderate  hardness  at  speeds 
and  obliquities  typical  of  current  ballistic  weapons.  The  first  stage 
of  this  investigation  is  restricted  to  predicting  trajectories  and  forces 
in  the  impact  of  spheres  against  thin  plates.  Spheres  are  chosen  because 
the  symmetry  of  the  sphere  simplifies  the  dynamics  of  the  interaction, 
thus  making  it  easier  to  evaluate  the  procedures  used  in  predicting  tra- 
jectories and  forces. 

The  predictive  procedure  consists  of  analytical  and  experimental 
parts.  The  analytical  procedure  predicts  trajectories  and  forces  from 
the  impact  speed,  obliquity,  radius  and  mass  of  the  sphere,  and  empirical 
constants  that  characterize  the  target.  Experimental  data  on  impacts  at 
normal  incidence  are  used  to  determine  the  empirical  constants  that  are 
needed  in  the  analytical  procedure. 


ANALYTICAL  PROCEDURE 


TYPICAL  IMPACT  PHENOMENA 

The  analytical  procedure  predicts  the  behavior  of  an  impact  system 
consisting  of  a very  hard  sphere  and  a significantly  softer  plate.  The 
prediction  is  made  by  representing  the  dynamics  of  the  real  system  by 
idealized  analytical  models  of  the  response  of  this  system  to  impact. 

In  order  to  explain  the  predictive  system,  the  principal  phenomena  of 
the  impact  of  a hard  sphere  against  the  softer,  comparatively  thin  plate 
will  be  described  in  this  section;  the  following  section  will  present 
approximations  and  assumptions  that  simplify  the  analytical  representa- 
tion of  these  phenomena. 

When  the  sphere  strikes  the  plate,  several  things  happen  simultan- 
eously: the  sphere  and  the  plate  begin  to  deform,  stresses  develop  in 
both  bodies  as  the  result  of  the  deformation,  and  both  bodies  are  accel- 
erated as  the  result  of  the  forces  on  the  surface  of  contact.  Interest 
is  here  restricted  to  impacts  that  leave  the  sphere  intact.  Such  impacts 
occur  if  the  sphere  is  very  hard  and  the  plate  is  much  more  ductile. 
Spheres  used  in  this  work  were  of  a hard  SAE  52100  steel,  and  the  plate 
was  a more  ductile  SAE  1010  steel  or  2024-T3  aluminum  alloy.  With  this 
system  the  spheres  were  recovered  intact  at  impact  speeds  up  to  1.2km/s. 
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Those  parts  of  the  craters  in  direct  contact  with  the  spheres  bore  the 
dimensions  and  curvature  of  the  spheres,  thus  demonstrating  that  the 
spheres  not  only  maintained  their  integrity  but  were  deformed  very  little 
during  impact,  while  the  plates  deformed  to  take  the  shape  of  the  sphere. 
Figure  1 is  an  example  of  crater  cross  sections  formed  by  essentially 
nondeforming  spheres. 

Observation  of  crater  shape  for  many  impacts  shows  that  the  sphere 
had  moved  along  curved  trajectories  for  all  impacts  except  those  at  nor- 
mal obliquity.  The  craters  and  high-speed  photographs  show  three  dis- 
tinct final  states  of  motion:  ricochet  from  the  front  surface  of  the 
plate,  embedment  V7ithin  the  plate,  and  perforation  of  the  plate.  In  the 
perforation  mode,  the  plate  fails  during  impact,  making  it  feasible  for 
the  sphere  to  pass  through  the  plate.  Failure  may  also  occur  in  either 
an  embedding  or  a ricocheting  impact,  but  it  does  not  influence  the  in- 
teraction enough  to  change  the  outcome.  The  observation  of  craters 
formed  by  oblique  impact  shows  that  all  such  impacts  begin  with  deflec- 
tions that  would  lead  to  ricochet.  Embedding  impacts  are  those  in  which 
the  projectile  speed  is  reduced  to  zero  during  the  interaction.  The  re- 
duction to  zero  velocity  could  occur  during  a ricochet  phase  of  the  im- 
pact or  after  a perforation  phase  has  begun. 

Plate  failure  at  near  normal  impact  obliquities  is  by  plugging,  i.e., 
by  separation  of  a cylindr ical-to-truncated  conical  plug  from  the  plate 
to  form  a symmetrical  hole  through  which  the  projectile  ultimately  passes 
(Figure  2).  At  larger  angles  of  obliquity  (roughly,  over  30  degrees), 
the  failure  process  loses  its  symmetry.  At  these  obliquities  perforation 
failures  are  different  for  aluminum  and  steel.  For  steel,  perforation 
occurs  by  the  formation  of  a fingerlike  piece  from  the  target  by  frac- 
tures that  are  cylindrical  at  one  end  and  with  two  fractures  on  each  side 
that  are  parallel  to  the  direction  of  motion  of  the  sphere  along  the  sur- 
face of  the  plate.  This  finger  bends  away  from  the  plate  and  frequently 
breaks  off.  When  it  breaks,  the  last  point  of  failure  is  at  the  hinge 
point  of  its  bending  (Figure  3a).  Failures  that  lead  to  perforation  in 
aluminum  are  markedly  different,  especially  for  thin  plates.  Failures 
in  thin  aluminum  plates  are  of  a petaling  type.  A central  fracture 
occurs  along  the  direction  of  travel  of  the  sphere,  and  additional  frac- 
tures travel  out  from  the  central  fracture  forming  petals  that  peel  back 
from  the  center  (Figure  3b) . 

RESTRICTIONS  ON  THE  IMPACT  SYSTEM  AND  APPROXIMATIONS 

The  analytical  procedure  used  to  represent  the  sphere-plate  system 
during  impact  is  based  on  a number  of  simplifying  restrictions  and  ap- 
proximations to  the  observed  behavior  of  the  system. 

First  it  is  assumed  that  the  analytical  procedure  can  be  projectile- 
oriented  and  can  represent  the  behavior  of  the  sphere-plate  system  ex- 
clusively in  terras  of  the  effect  on  the  motion  of  the  sphere.  Figure  4, 
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FIGURE  1.  Crater  Formed  by  a Sphere  That 
Remained  Intact.  Sphere:  SAE  52100  steel, 
6.35  mm  in  diameter;  plate:  SAE  1010  steel, 
1.47  mm  thick. 


which  illustrates  the  application  of  the  procedure  to  a ricochet  and  a 
perforation  trajectory,  helps  to  explain  the  role  of  restrictions  and 
approximations  in  the  analytical  procedure.  Thus,  both  for  a ricochet 
(Figure  4a)  and  for  a perforation  (Figure  4b),  the  trajectories,  which 
are  indicated  by  the  curves,  are  calculated  by  rigid-body  response  of 
the  sphere  to  forces  exerted  on  the  sphere  by  the  deforming  plate.  The 
deformations  of  the  plate  will  not  be  explicitly  determined,  and  it  is 
assumed  that  displacements  of  the  plate  material  are  such  as  to  accom- 
modate to  the  motion  and  shape  of  the  sphere.  In  this  process  of  de- 
formation the  plate  exerts  a force  on  the  sphere  that  depends  on  the 
motion  of  the  projectile  relative  to  the  plate  and  on  the  extent  of  con- 
tact between  the  sphere  and  the  plate.  This  point  of  view  is  permitted 
in  part  by  the  restriction  on  the  system  that  the  plate  and  its  support- 
ing structure  are  much  larger  and  more  massive  than  the  sphere,  so  that 
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FIGURE  2.  Cross  Section  of  a Crater  That 
Displays  Plugging-Mode  Failure.  Sphere: 

SAE  52100  steel,  6.35  nun  in  diameter;  plate: 
2024-T3  aluminum  alloy,  6.35  ram  thick. 


the  overall  motion  of  the  plate  is  not  significant  and  the  interaction 
of  the  plate  and  the  sphere  is  highly  localized  about  the  region  of  con- 
tact. Still  another  restriction  on  the  system  that  justifies  a com- 
pletely projectile-oriented  point  of  view  is  that  the  plate  material  is 
significantly  softer  than  the  projectile. 

The  simultaneous  and  interrelated  processes  of  deformation,  force 
development,  and  changes  in  motion  of  the  sphere  are  represented  in  sim- 
plified form  by  approximation  of  the  real  trajectory  of  the  sphere  by  a 
trajectory  composed  of  segments  having  constant  speed  and  curvature. 
These  segments  correspond  to  a finite  number  of  calculations  concerning 
location  of  the  sphere,  force  on  the  sphere,  and  change  in  its  state  of 
motion.  Figure  4 illustrates  ricochet  and  perforation  trajectories  that 
are  each  approximated  by  segments  of  constant  speed  and  curvature.  The 
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(a) 


(b) 


FIGURE  3.  Exit-Surface  Photographs  of  Two 
Craters  Formed  by  Impact  of  a 6.35-mm  SAE 
52100  Steel  Sphere  at  Obliquity,  Displaying 
(a)  Hinge-Mode  Failure  of  a 1.47-mmSAE  1010 
Steel  Plate,  and  (b)  Petaling-Mode  Failure 
of  a 1.27-mm  2024-T3  Aluminum  Alloy  Plate. 
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(a) 


(b) 
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FIGURE  4.  Hypothetical  Trajectories  Illustrating  the  Application 
of  the  Analytical  Procedure  to  (a)  Ricochet  and  (b)  Perforation. 
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basic  calculations  of  the  analytical  procedure  are  carried  out  at  the 
points  designated  by  the  heavy  dots. 

In  the  present  analytical  procedure  there  is  no  explicit  attempt  to 
do  detailed  calculations  of  the  strain  in  the  plate,  to  correlate  defor- 
mation directly  with  stress,  or  to  integrate  surface  stresses  and  trac- 
tions into  a resultant  force.  All  these  processes  are  implicit  in  the 
estimation  of  the  force  resisting  penetration  by  an  equation  of  the 
Poncelet  form  F = /l(a  + where  F is  the  resisting  force,  A is  the 

presented  area,  V is  the  impact  speed,  and  a and  h are  constants  that 
characterize  the  target.  The  dependence  of  the  force  on  deformation  is 
expressed  through  the  presented-area  factor,  and  the  dependence  on  speed 
is  given  by  the  quadratic  relation.  The  properties  of  the  system  are 
expressed  through  constants  a and  h.  The  resisting  function  given  above 
has  been  used  in  many  investigations  of  impact  at  normal  incidence . 

These  investigations  have  shown  that  this  function  is  capable  of 
predicting  the  outcome  of  impacts  for  a specific  system  with  reasonable 
accuracy.  Typically,  the  force  for  any  degree  of  embedment  is  estimated 
by  this  empirical  expression,  using  the  intersection  of  the  front  sur- 
face of  the  plate  and  the  surface  of  the  sphere  as  the  presented  area. 

In  order  to  determine  the  motion  of  the  sphere  at  oblique  impacts, 
some  rule  must  be  made  as  to  how  the  direction  of  the  force  is  to  be  de- 
termined from  the  geometry  of  the  contact  between  the  sphere  and  the 
plate.  The  initial  motion  of  the  sphere  is  in  the  plane  defined  by  the 
initial  line  of  flight  and  the  normal  to  the  plate.  The  force  vector 
that  governs  the  motion  after  contact  is  in  this  plane  because  of  the 
symmetry  of  the  contact  surface  about  this  plane.  The  rule  used  for  de- 
termining the  direction  of  the  force  within  this  plane  is  that  the  force 
is  normal  to  the  chord  connecting  the  upper  and  lower  extremities  of  the 
contact  surface  and  in  the  plane  of  motion.  This  chord  will  be  called 
the  contact  chovdy  and  the  plane  that  passes  through  the  line  of  chord 
and  is  perpendicular  to  the  plane  of  motion  will  be  called  the  contact 
plane.  The  presented  area  of  the  contact  surface  will  be  the  proj ection 
of  the  contact  surface  on  this  plane. 

The  force  on  the  sphere  is  considered  to  be  applied  to  the  center  of 
mass  of  the  sphere  in  accordance  with  the  principles  of  rigid-body  me- 
chanics. In  general,  there  must  be  a torque  applied  at  the  center  of 
mass  so  that  the  force  system  at  the  center  of  mass  is  equivalent  to  the 


^ F.  Helie.  Traite  de  Balestique  Expevimentale . Paris,  Dumaine,  1884. 
^ Ballistics  Institute.  Theory  of  Armor  Piercing  by  H.  Richter. 

St. -Louis,  France,  BI,  May  1946.  (Report  No.  3/46.) 

^ J.  Nishiwaki.  ^Resistance  to  Penetration  of  a Bullet  Through  an 
Aluminum  Plate,'*  J.  Phys.  Soc.  Jap.  y Vol.  6 (1951)  » P*  374. 
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force  system  applied  at  the  contact  surface.  In  the  present  procedure 
any  torque  of  this  sort  is  neglected  on  the  assumption  that  the  resulting 
rotation  and  angular  momentum  of  the  sphere  have  no  important  effect  on 
the  interaction.  The  force  at  the  center  of  mass  is  resolved  into  a 
force  component  in  the  direction  of  motion  and  a force  component  perpen- 
dicular to  the  direction  of  motion.  These  components  are  called  the  drag 
force  component  and  the  lift  force  component,  respectively,  using  ter- 
minology analogous  to  that  of  aerodynamics. 

At  some  point  of  penetration,  the  plate  fails.  For  any  given  calcu- 
lation this  can  be  assigned  a critical  penetration,  P*,  and  it  is  assumed 
in  the  analytical  procedure  that  for  any  system  there  is  a single  crit- 
ical penetration  for  failure  no  matter  what  the  obliquity  or  impact  ve- 
locity. It  is  further  assumed  that  the  only  important  effect  of  failure 
is  to  relieve  the  force  resisting  penetration  over  some  designated  part 
of  the  contact  surface.  The  extent  of  the  relieved  section  is  given  by 
an  angle  a*  that  is  measured  from  the  normal  to  the  plate  as  shown  in 
Figure  5.  The  effect  of  failure  on  the  trajectory  of  the  sphere  is  then 
automatically  taken  into  account  by  the  model  of  the  dynamics  of  the 
sphere . 

( 

DETAILS  OF  THE  ANALYTICAL  PROCEDURE 

This  section  summarizes  the  principal  steps  and  equations  that  are 
used  to  predict  the  results  of  impact.  This  prediction  procedure  re- 
quires as  input  data  the  mass,  m,  of  the  sphere;  its  radius,  its 

impact  speed,  VqI  the  obliquity,  0g;  the  thickness  of  the  plate,  T;  and 
the  constants,  a and  by  that  characterize  the  target  plate.  The  actual 
trajectory  of  the  sphere  is  approximated  by  a finite  number  of  segments, 
each  having  a constant  curvature  and  constant  speed  along  its  length. 

The  segment  represents  the  displacement  of  the  sphere  for  a known  time 
interval,  dt.  For  each  segment  the  main  steps  of  the  analytical  pro- 
cedure are  as  follows: 

1.  The  force  on  the  sphere  is  calculated  at  the  beginning  of  each 
segment  by  the  empirical  equation  and  with  direction  normal  to  the  con- 
tact plane. 

2.  The  force  is  resolved  into  a drag  component  that  is  in  the  dir- 
ection of  motion  and  a lift  component  that  is  in  the  direction  perpen- 
dicular to  the  motion  of  the  sphere. 

3.  The  change  in  speed  of  the  sphere  during  a time,  dty  is  calcu- 
lated from  the  drag  component  of  force,  and  the  curvature  of  the  path  is 
calculated  from  the  lift  component  of  force. 

4.  The  change  in  the  position  of  the  sphere  for  the  next  interval 
of  time,  dty  is  calculated  from  the  new  values  for  speed  and  curvature. 

5.  The  penetration  of  the  sphere  is  compared  with  the  critical  pen- 
etration for  failure  to  determine  whether  the  plate  has  failed  or  is 
intact. 
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FIGURE  5.  Diagram  of  the  Model  of  Failure  by  Plugging. 
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Force  on  the  Sphere  Before  Plate  Failure 

The  force  resisting  penetration  is  given  by  the  equation 

F = (a  + cos^  (1) 

where  F is  the  resisting  force,  A is  the  presented  area,  V is  the 
speed,  a and  b are  constants  that  characterize  the  target  plate,  and  ij; 
is  the  angle  between  the  sphere's  line  of  flight  and  the  normal  to  the 
contact  surface  (Figure  6). 

A general  expression  for  the  presented  area,  A^  is 


A = — ~(r  - sin  r ) cos  (()  + (r  - sin  F ) cos  (0  - ({))  (2) 


where  R is  the  radius  of  embedment, 
e 


= (2R  P - P-) 
e s 


1/2 


(3) 


the  angle  V is 


= 2ir  - 2 cos  ^[(R^  - P'>/(Rq  tan  6)] 


(4) 


and  when 


(F  - F)  / tan  0 > F , then  F = 2tt 
s e e 


The  angle  F is 

s 


-1 , 


= 2 cos  [(fl  - P)/(fl  sin  0)] 

o o o 


(5) 


and  when 


(F  - P)/sin  0 > F , then  F = 0 
s s s 


The  direction  of  the  force  is  specified  as  normal  to  the  contact 
plane;  thus,  the  force  is  at  an  angle,  ((),  with  respect  to  the  normal  to 
the  plate.  This  angle  is  defined  in  terms  of  the  obliquity,  0,  and  the 
penetration,  P,  as 

0 + sin  (P  - Rq)/Rq 

^ = (6) 

2 

and  this  angle  must  lie  between  zero  and  0, 
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VIEW  IN  THE 
DIRECTION  OF 
MOTION 


FIGURE  6,  Diagram  of  the  Forces  on  a Sphere.  Shaded  areas  are 
projected  contact  surfaces. 
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Model  of  Failure  of  the  Plate 


For  impacts  at  normal  incidence  the  failure  is  symmetrical,  and  a 
plug  is  immediately  expelled  from  the  plate.  Plugging  failure  is  put  in 
the  analytical  model  by  a penetration  criterion  that  determines  the  on- 
set of  failure.  The  expulsion  of  the  plug  is  represented  by  the  immedi- 
ate removal  of  resisting  force  over  the  region  of  the  plate  defined  by 
the  angle  a*.  Failure  of  the  plate  for  impacts  at  oblique  incidence  be- 
comes increasingly  asymmetrical  with  increasing  obliquity,  and  plate 
fractures  result  in  a petaling  formation  as  the  projectile  penetrates 
into  the  plate. 

The  model  of  failure  of  the  plate  for  impacts  at  obliquity  consists 
of  the  following  elements:  (1)  a critical  penetration,  P*,  for  the  onset 
of  failure;  (2)  an  angle,  a*,  to  define  the  initial  location  at  which 
failure  occurs;  and  (3)  a model  of  the  development  of  failure  as  the 
sphere  continues  its  motion.  There  are  two  options  here:  (1)  to  relieve 
resistance  immediately  as  for  normal  incidence,  and  (2)  to  maintain  con- 
tact over  a decreasing  contact  surface  so  as  to  model  the  peelback  of  a 
**hinge.  ** 

Elements  of  the  model  of  failure  are  illustrated  for  oblique  inci- 
dence in  Figure  7.  In  Figure  7a  the  sphere  has  reached  the  distance  P* 
at  which  failure  occurs.  The  location  of  the  failure  on  the  sphere  is 
given  by  the  angle,  a*,  which  is  measured  with  respect  to  the  normal  to 
the  exit  surface  of  the  plate,  so  that  the  specification  of  the  location 
of  failure  is  with  respect  to  the  plate  orientation.  In  the  first  option 
resistance  is  relieved  over  all  of  the  contact  surface  up  to  the  angle 
a*.  In  Figure  7B  the  sphere  has  continued  its  trajectory,  and  the  dif- 
ference between  angles  U2  and  defines  the  opening  of  the  failure. 
These  angles  are  computed  on  the  premise  that  the  lengths  of  the  plate 
elements  in  contact  with  the  sphere  remain  invariant  above  and  below  the 
failure  as  the  sphere  moves  and  the  failure  opens  up.  This  is  done  by 
keeping  the  angle  Q constant  for  the  upper  segment  and  by  reducing  the 
lower  segment  by  a length  increment,  dS\  for  an  increment  of  motion,  d5, 
and  change  of  obliquity,  d0,  that  is  given  by 

dS^  = dS  ^ RdQ  (7) 

o 

where  R^dQ  is  the  swing  of  the  radius  of  the  sphere  through  the  incre- 
ment of  obliquity,  dd. 


Force  on  the  Sphere  After  Plate  Failure 


After  the  plate  has  failed,  the  calculation  of  force  has  to  account 
for  relief  of  the  force  over  part  of  the  plate  due  either  to  plugging  or 
petaling.  For  symmetrical  failure  by  plugging,  the  removed  part  of  the 
plate  is  defined  by  the  intersection  of  the  spherical  surface  by  a plane 
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FIGURE  7.  Diagram  of  the  Hinge  Mode  of  Failure. 
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that  is  parallel  to  the  front  surface  of  the  plate  and  displaced  from  it 
by  a distance,  /z,  that  is  related  to  the  penetration  and  angle  a by  the 
equation 

h = R cos  a + P - i?  (8) 

s s 

Procedures  for  determining  the  presented  area  and  the  inclination  of  the 
contact  surface  are  identical  to  those  used  for  the  calculation  of  force 
before  the  failure  of  the  plate  (Equations  1 through  6) , except  that  the 
plane  of  reference  is  displayed  by  the  distance  h (Figure  8).  The  equa- 
tions for  the  presented  area  are 

i?2 

i4'=-^  (r^  - sin  rpcos  (j)'  +—  (r^  - sin  rpcos  (0  - (}>')  (9) 

R^=  [2R^(P  - h)  - iP  - ?z)2]  1^2 

where 

= 2tt  - 2 cos"^  [(i?g  - P + h)/{R^  tan  0)]  (11) 

r'  = 2 cos"^  [(P  - P + h)/(R  sin  0)]  (12) 

s s s 

and  for  inclination  of  the  contact  plane, 

0 + sin"\p  - h - R )/R 

<!>' ^ ^ (13) 

2 

The  force  calculated  for  this  area  is 


F'  = A' (a  + bV^  cos2  ip)  (14) 

Let  P*  be  the  same  as  F of  Equation  1;  then  the  net  force,  F,  acting  on 
the  sphere  after  plugging  failure  is  the  force  F*  calculated  without 
plate  failure  minus  the  force  F',  the  force  relieved  by  the  plugging 
failure;  i. e. , 

F = (15) 

The  analytical  model  of  the  hinge  mode  of  failure  uses  the  relieving 
effect  of  the  plugging  model  but  in  addition  includes  the  resistance  due 
to  the  finger  that  is  formed  by  the  unsymmetrical  failure  characteristic 
of  oblique  impact.  This  resistance  is  calculated  as 
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FIGURE  8.  Diagram  of  the  Forces  on  a Sphere  After  Plate  Failure. 
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F"  = A^ia  + bV^  cos2  (16) 

where 

A = n[i?  sin  (8/2)]^  (17) 

o S 

and  3 is  the  angle  between  the  upper  edge  of  the  hinge  and  the  last  point 
of  contact  with  the  sphere,  as  shown  in  Figure  9.  This  angle  begins  with 
the  value  g = a + 90°  - 0 if  this  value  is  less  than  2a;  if  not,  then 
3 = 2a  until  0 is  greater  than  0^.  ^ For  each  segment  of  the  path  there- 
after, the  value  of  3 is  corrected  for  the  displacement  of  the  segment 
by 


3'  = 3 - do  - (dS/R)  (18) 

o 

The  angle  ' is  calculated  by 

= 90°  + 0 - 0^  - (3/2)  (19) 

if  0 - 0^  is  negative,  and  by 

= 90°  - (3/2)  (20) 

if  the  difference  is  positive. 

The  net  force  under  conditions  of  hinge  failure  is 

F = F*  - F'  + (21) 


Dynamics  and  Kinematics  of  the  System 

The  dynamics  of  the  impact  of  the  sphere  against  the  plate  are 
modeled  by  calculating  the  force  on  the  sphere  by  the  method  just  de- 
scribed and  then  calculating  the  rigid-body  response  of  the  sphere  to 
this  force.  It  is  convenient  to  resolve  the  force  into  a component  that 
is  in  the  current  direction  of  motion  and  a component  that  is  perpendic- 
ular to  the  current  direction  of  motion.  Changes  in  the  speed  of  the 
sphere  are  calculated  from  the  component  that  is  in  the  direction  of  mo- 
tion, and  the  curvature  of  the  trajectory  is  calculated  from  the  compon- 
ent perpendicular  to  the  direction  of  motion. 

The  force  calculated  by  Equations  1,  15,  or  21  is  at  an  angle, 
with  respect  to  the  direction  of  motion  of  the  sphere.  This  angle  can 
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FIGURE  9 


Diagram  of  the  Parameters  of  Hinge  Failure. 
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be  determined  from  the  obliquity,  0,  and  the  inclination  of  the  contact 
plane,  (J): 


Ij;  = 0 - 

<)) 

(22) 

The  drag  and  lift  components  of  the  force  are  then 

II 

cos 

(23) 

II 

sin 

(24) 

The  changes  in  velocity  and  the  curvature  of  the  trajectory  correspond- 
ing to  these  components  of  force  are 

dV  = F^dt/m  (25) 

and 

r = mV'^lF^  (26) 

L 


where  dt  is  the  increment  of  time  over  which  the  calculation  is  carried 
out,  and  r is  the  radius  of  curvature  of  the  trajectory, 

These  quantities,  dV  and  r,  determine  the  segments  of  the  trajector- 
ies. The  obliquity  is  changed  as 

de  = Vdtiv  (27) 

so  that  the  new  obliquity  is  0', 

0 ' = 0 + J0  (28) 


and  the  corresponding  increment  of  penetration  is 


dP  = V cos  0dt 


P'  = P + dP 


(29) 


and  the  corresponding  increment  of  lateral  motion  is 

dY  ^ V sin  0dt  y'  = y + dy  (30) 

These  steps  are  illustrated  in  Figure  10,  where  it  should  be  noted  that 
segment  size  is  exaggerated  to  facilitate  the  explanation. 
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FIGURE  10.  Diagram  of  the  Kinematics  of  Trajectories, 
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Prediction  of  Crater  Shape 

The  sphere  sweeps  out  a tube  as  it  moves  along  its  trajectory.  The 
crater  in  the  plate  is  the  part  of  the  volume  of  the  tube  extending  be- 
yond the  front  surface  of  the  plate.  The  shape  of  the  crater  is  speci- 
fied completely  by  its  cross  section  in  the  plane  of  motion  of  the 
sphere;  thus,  crater  shapes  can  be  calculated  from  trajectories  as  the 
projection  on  the  X^Y  plane  of  the  tube  intersection  that  is  swept  out 
by  the  sphere,  and  the  half  space  defined  by  the  plane  of  the  front  sur- 
face of  the  plate.  The  projection  of  the  tube  is  two  curves.  The  curve 
corresponding  to  the  side  of  the  projection  defining  the  crater  shape  is 
determined  from  the  trajectory  as  the  focus  of  points  given  by 


X = X ■¥  R sin0 
s 

(31) 

y = y ~ i?  cos  0 
^ s 

(32) 

where  x and  y are  points  on  the  projection  of  the  tube,  X and  Y are 
points  on  the  trajectory,  and  R and  0 are,  respectively,  the  radius  of 
the  sphere  and  the  obliquity  of  the  motion  of  the  sphere  at  a given  point 
on  its  trajectory.  The  intersection  of  this  curve  with  the  half  space 
defined  by  the  front  surface  of  the  plate  is  the  curve  that  describes 
the  crater  shape. 

It  is  a minor  convenience  to  define  x and  y for  each  calculation  in 
such  a way  that  x and  y are  both  zero  at  the  entry  intersection  of  the 
crater  shape  curve  with  the  front  surface  of  the  plate.  Trajectories 
are  arbitrarily  defined  to  have  initial  values 


X = -R 

s 

(33) 

y = 0 

(34) 

The  corresponding  values  of  x and  y are 

X = -R  + i?  sin  0 

s s 

y - cos  0 (36) 

In  general,  y will  not  be  zero  when  x is  zero.  A set  of  coordinates  x', 
y'  such  that  the  entry  intersection  at  the  front  surface  of  the  plate  is 
given  by  x'  = 0 and  y ' = 0 is  determined  in  our  finite  step  procedure  by 
the  simple  transformation 
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x'  = X (37) 

y’  = y + (38) 

where  y is  determined  from  y,,  the  y value  corresponding  to  the 

value  o?  X that  is  positive,  By  the  equation 

^o  ° ^1  ■ ^1  ® 

EXPERIMENTAL  DETERMINATION  OF  EMPIRICAL  CONSTANTS 

The  analytical  model  requires  the  mass  and  radius  of  the  sphere  and 
the  thickness  of  the  plate  as  basic  parameters  to  define  the  size,  con- 
figuration, and  relevant  inertial  properties  of  the  system.  In  addition, 
the  model  requires  four  empirical  parameters  that  account  for  mechanical 
properties  of  the  plate.  The  constants  a and  b of  the  force  equations 
represent  the  capacity  of  the  plate  to  exert  force  on  the  sphere  as  the 
plate  deforms.  The  parameter  P*  determines  the  beginning  of  failure, 
and  the  parameter  a*  its  location.  These  four  parameters  must  be  deter- 
mined either  by  separate  theoretical  considerations  or  empirically.  In 
the  present  work  these  are  determined  empirically. 

The  parameters  P*  and  a*  can  be  estimated  directly  from  observations 
on  impact  craters.  A series  of  observations  on  impacts  with  increasing 
impact  speed  shows  that  failure  occurs  near  a constant  value  of  penetra- 
tion. When  failure  does  occur,  it  is  possible  to  fix  the  location  on  the 
surface  of  the  sphere  and  estimate  the  angle  a*.  Such  observations  on 
craters  are  basically  consistent  with  our  assumption  that  failure  can  be 
described  by  assigning  constant  values  of  P*  and  a*.  However,  there  are 
observable  variations  in  these  parameters  as  the  speed  and  obliquity  are 
varied,  so  that  constancy  of  these  parameters  as  a description  of  fail- 
ure is,  at  best,  a good  approximation.  For  example,  there  are  indica- 
tions that  the  angle  a*  tends  to  decrease  with  increasing  impact  speed, 
and  although  the  penetration  P*  appears  to  be  quite  constant  for  alum- 
inum, it  varies  considerably  for  mild  steel. 

Constants  P*  and  a*  are  selected  to  match  observations  of  crater 
formation  at  normal  incidence;  then  the  constants  a and  i>,  which  char- 
acterize the  target,  are  derived  from  experimental  data  in  a two-step 
procedure.  A statistical  procedure  operates  on  experimental  data  on  the 
residual  speed  as  a function  of  the  impact  speed  at  normal  incidence  to 
determine  the  best  values  of  the  constants  R and  V in  the  DRI  equation^ 


7 = P /7^  - 72  (40) 

X s o 


^ R.  Recht  and  T.  W.  Ipson.  "Ballistic  Perforation  Dynamics,"  J, 
Appt.  Meah. , Vol.  30,  Series  E,  No.  3 (September  1963). 
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An  iterative  procedure  determines,  by  trial  and  correction,  the  em- 
pirical constants  a and  h that  give  the  best  agreement  between  calcula- 
tions done  at  a predetermined  number  of  impact  speeds  and  the  predictions 
of  the  above  empirical  formula. 

For  a given  system,  data  on  residual  speed  as  a function  of  impact 
speed  for  impacts  at  normal  incidence  is  fit  to  Equation  40  by  a least- 
squares  fit  of  the  data  to  the  linear  expression 

Y = MX  ^ B (41) 

where  Y = 7^,  X = 7^,  M = y and  B = i?^7^.  Trial  values  of  a and  b are 
then  selected  and  used  in  the  analytical  model  for  impacts  at  normal  in- 
cidence to  obtain  predictions  of  residual  speeds  that  are  compared  with 
the  empirical  equation.  An  iterative  procedure  then  introduces  small 
changes  in  a and  b systematically,  so  that  differences  between  the  pre- 
dicted results  and  experimental  results  are  minimized. 


trial  values  are 

a = mV^/luR^ 
o s 

(42) 

b = wR^/l'nR^ 
s 

(43) 

These  values  are  used  in  the  analytical  procedure  for  normal  incidence 
to  determine  the  residual  speed  for  several  impact  speeds.  These  pairs 
of  residual  speeds  and  impact  speeds  are  then  used  to  determine  a pair 
of  values  of  R and  7^  by  the  same  curve-fitting  procedure  that  was  used 
on  the  empirical  data.  Constant  a is  varied  by  a predetermined  value, 
dUy  and  another  pair  of  values  for  R and  7q  is  calculated.  Then  b is 
varied  by  a predetermined  amount,  dby  and  a new  pair  of  values  of  R and 
7q  is  calculated.  From  the  information  available  at  this  point,  four 
partial  derivatives  are  approximated  by  finite  differences  as 


dV  7 (a  + da,b)  - 7 (a,b) 

o O O '’ 

da  da 


(44) 


dR  _ R(a  + da^b)  - R(a^b) 
da  ~ da 


(45) 


97  7 (a^b  + db)  - 7 (a^b) 

o _ o o ^ ^ 

db  ~ db 


(46) 
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9i?  _ R(a,b  + db)  - R(atb) 
db  " db 


These  partial  derivatives  and  the  differences  dV^  = (experimental) 

- V (cijb)  and  dR  = R (experimental)  - R (a^b)  are  used  to  determine  new 
values  of  dsc  and  db  by  solving  the  equations 


3K  91/ 


and 


dR  = 


9a 


da 


j.  Ah 


(48) 


(49) 


and  the  process  is  repeated  until  the  differences  are  below  some  preset 
value. 

In  summary,  this  model  of  impact  dynamics  requires  basic  parameters 
that  define  the  configuration,  size,  and  inertia  of  the  system.  In 
addition,  this  model  uses  four  empirical  constants  that  define  the  ma- 
terial response  of  the  plate  to  impact.  Two  empirical  constants,  P*  and 
a*,  are  determined  from  observations  on  the  depth  of  penetration  required 
for  the  failure  of  the  plate  and  on  location  of  the  failure.  The  other 
two  constants,  a and  Z?,  are  determined  by  a best  fit  to  data  on  the  re- 
sidual speed  as  a function  of  impact  speed  for  normal  incidence.  Cal- 
culations of  residual  speeds,  exit  obliquities,  and  trajectories  of  the 
sphere  through  the  plates  are  compared  directly  with  experimental  data. 
The  degree  of  agreement  is  taken  as  a measure  of  success  in  modeling  the 
dynamics  of  the  impact.  The  agreement  is  very  good  at  normal  incidence 
where  it  is  essentially  a fit  to  that  data,  and  agreement  merely  shows 
that  it  is  possible  to  fit  experimental  data  by  means  of  the  given  pro- 
cedure. The  degree  of  agreement  at  other  obliquities  indicates  success 
in  extending  the  model  of  impact  dynamics  to  the  more  complicated  condi- 
tions of  oblique  impact. 


COMPARISON  OF  ANALYTICAL  PREDICTIONS  AND  EXPERIMENTAL  DATA 


The  experimental  part  of  this  work  consisted  of  experimental  firings 
of  spheres  against  steel  and  aluminum  alloy  plates.  The  firings  were 
covered  by  high-speed  photographic  instrumentation,  electronic  measure- 
ment of  impact  speed,  and  measurement  of  craters  produced  in  the  target 
plates.  The  final  form  of  the  data  is  the  correlation  of  the  exit  speed 
and  obliquity  of  the  sphere  with  its  initial  speed  and  obliquity,  and 
measurement  of  the  craters  that  were  formed  at  particular  impact  speeds 
and  obliquities. 
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Chosen  for  this  experimental  work  were  6.35-min  spheres  of  SAE  52100 
alloy  steel.  The  target  plates  were  6061-T6  and  2024-T3  aluminum  alloy, 
and  a mild  steel  plate  (SAE  1010).  The  aluminum  alloy  plates  varied  from 
1.27  to  9.53  mm  in  thickness.  The  steel  plates  were  all  1.47  mm  thick. 
Firings  were  performed  with  the  plates  at  normal  incidence  to  the  path 
of  the  projectile  in  order  to  collect  sufficient  data  to  determine  the 
constants  needed  for  the  expression  for  the  force  on  the  sphere.  Table  1 
lists  the  values  of  a,  P*,  and  a*  that  were  determined  from  the  data 
at  normal  incidence.  Additional  firings  at  oblique  incidence  were  per- 
formed at  30,  45,  and  60  degrees  as  a rule,  but  a few  data  were  collected 
for  20,  50,  and  70  degrees. 


TABLE  1.  Constants  Determined  From  Data  at  Normal  Incidence. 


Target  plate 

Target  plate 
thickness , 
mm 

Constants 

material 

pt  mm 

a*  deg 

a,  N/m2 

fc,  kg/m^ 

2024-T3  aluminum 
alloy 

1.27 

2.2 

45 

0.272  X 10^ 

0.129  X lO'* 

3.18 

3.0 

45 

0.101  X 10^° 

0.126  X lO** 

6.35 

6.4 

45 

0.933  X 10^ 

0.974  X 103 

9.52 

7.5 

45 

0.112  X 10^0 

0.580  X 103 

SAE  1010  steel. . . 

1.47 

4.6 

45 

0.380  X 10^ 

0.148  X lO'* 

A given  sphere  and  plate  constitute  a system.  Correlation  of  the 
final  state  of  motion  of  the  sphere  to  the  initial  state  of  motion  is  a 
measure  of  the  impact  behavior  of  the  system.  The  accuracy  of  the  pre- 
diction of  this  correlation  is  taken  as  an  index  of  the  accuracy  of 
modeling  of  the  system  during  impact.  Impact  systems  of  the  same  ma- 
terial but  different  plate  thicknesses  constitute  a family  of  systems. 
The  case  of  the  steel  sphere  and  aluminum  plates  has  been  investigated 
as  a family,  and  it  is  possible  to  come  to  conclusions  about  the  family, 
especially  as  to  the  meaning  of  the  constants  for  the  force  expression 
as  a function  of  plate  thickness. 


STEEL  SPHERES  IMPACTING  AGAINST  ALUMINUM  PLATES 


Figure  11  shows  experimental  data  compared  to  the  analytical  repre- 
sentation for  firings  of  6.35-mm  steel  spheres  against  9.53-mm  aluminum 
plates  at  normal  incidence.  The  data  points  are  the  results  of  experi- 
mental firings.  The  line  represents  the  analytical  predictions  based 
on  the  model  described  above  and  the  values  of  constants  a and  h in  the 
expression  for  the  contact  force  that  give  the  best  fit  to  the  data. 
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FIGURE  11.  Data  and  Predictions  of  Exit  Speeds  for  6.35-mm 
Spheres  Fired  Against  6061  Aluminum  Alloy  Targets  of  9.53- 
mm  Thickness  at  Zero  Degree  Obliquity. 


This  comparison  is  only  an  indication  of  the  consistency  of  the  analyt- 
ical procedure  with  the  experimental  data.  It  demonstrates  the  feasibil- 
ity of  finding  constants  for  the  empirical  expression  for  the  contact 
force  in  the  analytical  model,  and  results  in  predictions  consistent 
with  data  used  to  derive  the  constants. 

Agreement  between  predictions  at  obliquity  and  experimental  data  is 
considerably  more  significant  than  the  agreement  at  normal  incidence  be- 
cause there  are  qualitative  features  of  the  data  predicted  by  the  model 
that  are  not  obvious  from  the  data  at  normal  incidence.  The  system  with 
the  largest  plate  thickness  is  chosen  first  because  this  system  displays 
these  qualitative  features  in  a particularly  simple  form. 
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Thick-Plate  Ricochet 


Figures  12  through  17  show  a comparison  of  experimental  data  and 
analytical  predictions  for  firings  of  6.35-mm  spheres  against  6061-T6 
aluminum  alloy  targets  at  20,  45,  and  60  degrees  obliquity.  For  each 
obliquity  there  is  a range  of  impact  speeds  in  which  the  sphere  rico- 
chets, a range  in  which  it  embeds,  and  a range  in  which  it  perforates 
(see  Figure  12).  There  is  reasonable  agreement  between  the  experimental 
data  and  the  predictions  as  to  the  extent  of  these  ranges.  The  curves 
of  residual  speed  in  the  ricochet  range  are  roughly  sinusoidal.  The 
corresponding  ricochet  obliquities  increase  monotonically  from  a value 
above  90  to  180  degrees  and  then  remain  constant.  This  is  what  the  model 
predicts,  and  the  experimental  data  in  general  agree  with  these  predic- 
tions. In  the  embedment  range  the  residual  speed  is  zero,  but  the  ana- 
lytical predictions  assign  an  obliquity  to  every  impact  speed  in  the 
embedment  range.  The  obliquity  assigned  is  the  last  obliquity  angle  the 
projectile  has  as  it  comes  to  rest.  In  terms  of  the  trajectory  of  the 
sphere,  it  is  the  terminal  direction  of  the  trajectory.  Examples  of 
these  predictions  are  shown  in  the  embedment  zones  of  Figures  13,  15,  and 
17.  In  the  perforation  range  the  residual  speed  curves  are  of  gradually 
decreasing  curvature,  and  the  corresponding  curves  for  the  obliquity 
after  perforation  rise  rapidly  and  approach  the  incident  obliquity. 

Figures  18  through  24  illustrate  the  same  kind  of  comparison  of  pre- 
dictions and  data  for  a system  consisting  of  a 6.35-mm  sphere  and  2024-T3 
aluminum  alloy  plates  of  6.35-mm  thickness.  The  boundaries  of  the  embed- 
ment, ricochet,  and  perforation  ranges  are  changed,  but  the  qualitative 
descriptions  of  the  correlations  in  each  range  remain  the  same.  The 
system  with  the  9.53-mm  6061-T6  aluminum  alloy  plates  and  the  system  with 
the  6.35-mm  2024-T3  aluminum  alloy  plates  both  have  the  important  char- 
acteristic that  the  plate  fails  at  deep  penetrations. 

Consider  the  trajectories  of  a projectile  in  the  ricochet  and  embed- 
ment range.  As  impact  speed  increases,  the  maximum  penetration  increases. 
This  is  predicted  by  the  model  and  observed  in  experimental  firings.  If 
the  plate  is  sufficiently  thick,  the  sphere  will  never  achieve  the  pene- 
tration required  for  the  failure  of  the  plate  on  a trajectory  that  is  of 
the  ricochet  type.  Instead,  failure  will  occur  on  a trajectory  that  ends 
in  embedment  even  though  the  plate  has  failed  at  some  point  along  the 
trajectory. 

These  two  systems  exhibit  what  may  be  called  "thick-plate  ricochet 
behavior."  The  correlation  of  exit  speed  with  impact  speed  is  smooth 
and  S3nmnetrical.  The  correlation  of  exit  obliquity  with  impact  speed  is 
a smooth  increase  to  a flat  top.  These  are  the  characteristics  of  rico- 
chet with  no  tear  or  fracture  failures  that  go  completely  through  the 
plate. 
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FIGURE  12.  Data  and  Predictions  of  Exit  Speeds  for  6.35-tnin  Spheres 
Fired  Against  6061  Aluminum  Alloy  Targets  of  9.53-iran  Thickness  at  20 
Degrees  Obliquity. 
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FIGURE  13.  Data  and  Predictions  of  Exit  Obliquities  for 
6.35-mm  Spheres  Fired  Against  6061  Aluminum  Alloy  Targets 
of  9.53-mm  Thickness  at  20  Degrees  Obliquity. 
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FIGURE  14.  Data  and  Predictions  of  Exit  Speeds  for  6 . 35-mm  Spheres 
Fired  Against  6061  Aluminum  Alloy  Targets  of  9.53-mm  Thickness  at  45 
Degrees  Obliquity. 
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FIGURE  15.  Data  and  Predictions  of  Exit  Obliquities 
Spheres  Fired  Against  6061  Aluminum  Alloy  Targets  of 
ness  at  45  Degrees  Obliquity. 


for  6.35-mm 
9 . 53-rara  Thick- 
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FIGURE  16.  Data  and  Predictions  of  Exit  Speeds  for  6.35-mm  Spheres 
Fired  Against  6061  Aluminum  Alloy  Targets  of  9.53-mm  Thickness  at 
60  Degrees  Obliquity. 
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FIGURE  17.  Data  and  Predictions  of  Exit  Obliquities  for  6.35-mm 
Spheres  Fired  Against  6061  Aluminum  Alloy  Targets  of  9 . 53-mm  Thick- 
ness at  60  Degrees  Obliquity. 
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FIGURE  18.  Data  and  Predictions  of  Exit  Speeds  for  6.35-min  Spheres 
Fired  Against  2024  Aluminum  Alloy  Targets  of  6.35-mm  Thickness  at 
Zero  Degree  Obliquity. 
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FIGURE  19.  Data  and  Predictions  of  Exit  Speeds  for  6 . 35-mm  Spheres 
Fired  Against  2024  Aluminum  Alloy  Targets  of  6.35-mm  Thickness  at 
30  Degrees  Obliquity. 
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FIGURE  20.  Data  and  Predictions  of  Exit  Obliquities  for  6.35-mm 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  6 . 35-mm  Thick- 
ness at  30  Degrees  Obliquity. 
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FIGURE  21.  Data  and  Predictions  of  Exit  Speeds  for  6.35-nnn  Spheres 
Fired  Against  2024  Aluminum  Alloy  Targets  of  6.35-mm  Thickness  at 
45  Degrees  Obliquity. 
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FIGURE  22.  Data  and  Predictions  of  Exit  Obliquities  for  6.35-mm 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  6 . 35-mm  Thick- 
ness at  45  Degrees  Obliquity. 
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FIGURE  23.  Data  and  Predictions  of  Exit  Speeds  for  6.35-n\m  Spheres 
Fired  Against  2024  Aluminum  Alloy  Targets  of  6.35-mm  Thickness  at 
60  Degrees  Obliquity. 
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FIGURE  24.  Data  and  Predictions  of  Exit  Obliquities  for  6. 35-ram 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  6 . 35-mm  Thick- 
ness at  60  Degrees  Obliquity. 


Thin-Plate  Ricochet 


Figures  25  through  34  show  behavior  of  a system  consisting  of  a 6.35- 
iran  sphere  and  either  a 3.18-  or  a 1.27-mm  2024-T3  aluminum  alloy  plate. 
These  correlations  are  similar  to  those  for  the  thicker  plates.  The  most 
conspicuous  differences  are  that  correlations  of  exit  speed  with  impact 
speed  are  not  symmetrically  shaped,  and  correlations  of  exit  obliquity 
with  impact  speed  have  pronounced  spikes  near  the  embedment  range.  These 
differences  are  caused  by  ultimate  failure  of  the  plate  in  the  ricochet 
range.  Failure  occurs  along  a trajectory  that  ends  in  ricochet  without 
failure.  Failure  influences  the  trajectory,  but  not  enough  to  result  in 
either  embedment  or  perforation.  The  effect  is  rather  to  ricochet,  but 
at  a lower  exit  speed  and  a higher  exit  obliquity.  There  is  a small 
range  of  impact  speeds  for  which  this  modified  form  of  ricochet  occurs. 
For  speeds  above  this  range  the  trajectories  end  in  embedment. 

In  a given  family  of  impact  systems,  this  change  in  behavior  is  cor- 
related with  thickness,  and  below  some  critical  thickness  the  systems  of 
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FIGURE  25.  Data  and  Predictions  of  Exit 
Speeds  for  6.35-mm  Spheres  Fired  Against 
2024  Aluminum  Alloy  Targets  of  3.18-mm 
Thickness  at  Zero  Degree  Obliquity. 


a family  will  always  show  this  kind  of  behavior.  The  controlling  physi- 
cal phenomena  are  the  processes  of  failure  in  the  plate  that  take  place 
at  lower  velocities  as  the  plate  thickness  is  decreased.  The  dashed 
lines  in  Figures  30  and  31  are  analytical  predictions  using  an  extremely 
large  value  of  the  critical  penetration.  This  shows  the  thick-plate  be- 
havior of  the  system  compared  with  the  actual  thin-plate  behavior  of  the 
system. 

The  comparison  between  experimental  data  and  analytical  prediction 
leads  to  the  conclusion  that  the  analytical  model  is  sufficiently  accur- 
ate to  fit  data  at  normal  incidence  by  means  of  its  empirical  constants. 
In  addition,  it  can  predict  qualitative  features  of  the  results  of  impact 
at  oblique  incidence  that  are  not  evident  from  observations  at  normal 
incidence.  Such  features  as  the  spikes  in  the  correlations  of  exit  ob- 
liquities with  impact  speed  and  the  rapid  drops  in  the  correlations  of 
exit  speed  with  impact  speed  are  unforeseen  consequences  of  the  model  of 
impact  dynamics.  The  existence  of  these  features  is  predicted  and  found 
to  be  in  reasonable  quantitative  agreement  with  experimental  observa- 
tions . 
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FIGURE  26.  Data  and  Predictions  of  Exit  Speeds  for  6.35-inin 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  3.18-mm 
Thickness  at  30  Degrees  Obliquity. 
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FIGURE  27.  Data  and  Predictions  of  Exit  Obliquities  for 
6.35-mm  Spheres  Fired  Against  2024  Aluminum  Alloy  Targets 
of  3.18-mm  Thickness  at  30  Degrees  Obliquity. 
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FIGURE  28.  Data  and  Predictions  of  Exit  Speeds  for  6.35-nim 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  3.1&-mm 
Thickness  at  45  Degrees  Obliquity. 
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FIGURE  29.  Data  and  Predictions  of  Exit  Obliquities  for 
6.35-iran  Spheres  Fired  Against  2024  Aluminum  Alloy  Targets 
of  3.18-mm  Thickness  at  45  Degrees  Obliquity. 
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FIGURE  30.  Data  and  Predictions  of  Exit  Speeds  for  6.35-mm 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  3.18-mm 
Thickness  at  60  Degrees  Obliquity. 
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FIGURE  31.  Data  and  Predictions  of  Exit  Obliquities  for  6.35- 
nun  Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  3.18- 
mm  Thickness  at  60  Degrees  Obliquity. 
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FIGURE  32.  Data  and  Predictions  of  Exit  Speeds  for  6.35-mm 
Spheres  Fired  Against  2024  Aluminum  Alloy  Targets  of  1.27“ 
mm  Thickness  at  Zero  Degree  Obliquity. 
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FIGURE  33.  Data  and  Predictions  of  Exit 
Speeds  for  6.35-mm  Spheres  Fired  Against 
2024  Aluminum  Alloy  Targets  of  1.27-mm 
Thickness  at  60  Degrees  Obliquity. 
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FIGURE  34.  Data  and  Predictions  of  Exit 
Obliquities  for  6.35-min  Spheres  Fired 
Against  2024  Aluminum  Alloy  Targets  of 
1.27-iran  Thickness  at  60  Degrees  Obliquity. 
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Crater-Shape  Predictions 

Correlations  between  initial  and  final  states  of  motion  are  useful 
representations  of  system  behavior,  and  the  accuracy  with  which  a model 
predicts  these  correlations  is  one  measure  of  the  accuracy  of  the  model 
for  that  system.  Prediction  of  crater  shapes  for  given  impact  speeds  is 
an  even  more  sensitive  measure  of  the  accuracy  of  a model. 

Figure  35  compares  actual  crater  shapes  with  predictions  from  the 
analytical  model.  The  solid  lines  represent  tracings  of  craters  taken 
from  a plate  that  has  been  cut  through  the  crater  in  the  direction  of 
travel  of  the  sphere.  Dashed  lines  represent  predictions  derived  from 
the  trajectory  of  the  sphere  by  applying  Equations  31,  32,  37,  and  38. 

The  impact  system  is  a 6.35-mm  sphere  impacting  a 9.53-mm  aluminum  alloy 
plate  at  60  degrees  obliquity.  The  craters  correspond  to  points  in  the 
ricochet  and  embedment  ranges  of  Figures  16  and  17,  the  correlation  of 
initial  with  final  states  of  motion  for  this  impact  system. 

The  predictions  and  tracings  have  roughly  similar  shapes,  but  a sig- 
nificant difference  in  major  dimensions.  Actual  craters  are  shorter  and 
shallower  than  the  predicted  craters  in  the  lower  part  of  the  ricochet 
range.  With  increasing  impact  speed,  but  remaining  in  the  ricochet 
range,  the  actual  crater  length  exceeds  the  predicted  crater  length, 
and  the  actual  crater  shape  is  more  curved  toward  the  exit  end.  In  the 
embedment  range,  actual  and  predicted  shapes  are  in  reasonable  agreement. 
Perforation  craters  offer  the  least  evidence  of  the  dynamics  of  the  im- 
pact process  because  large  sections  of  the  craters  are  torn  away. 

Figure  36  compares  tracings  of  crater  shapes  with  predictions  for 
the  system  consisting  of  a 6.  35 -mm  sphere  and  a 6.35-’iran  aluminum  alloy 
plate.  These  comparisons  have  the  same  trends  as  for  the  9.53-mm  6061- 
T6  plates,  but  the  difference  between  predicted  and  observed  crater 
shapes  is  greater.  This  is  consistent  with  the  observation  that  predic- 
tions of  the  final  states  of  motion  are  poorer  for  this  plate  as  well. 

Figure  37  compares  tracings  of  crater  shapes  with  predictions  for  a 
system  consisting  of  a 6.35-mm  sphere  and  a 3.18-nnn  2024-T3  aluminum 
alloy  plate.  Comparisons  for  this  system  are  considerably  closer  than 
for  the  6.35-mm  plate  and  show  the  same  trends — an  overprediction  at 
lower  speeds,  and  improved  agreement  with  increased  impact  speed. 

Steel  Spheres  Impacting  Against  Mild  Steel  Plates 

It  is  not  practical  to  conduct  a study  with  steel  plates  comparable 
to  the  study  with  aluminum  plates.  The  difficulty  is  that  the  greater 
density  and  strength  of  steel  results  in  breakup  of  the  sphere  at  mod- 
erate plate  thicknesses.  A series  of  experiments  has  been  conducted 
with  1.47-mm-thick  mild  steel,  and  these  data  have  revealed  still  other 
complexities  in  the  behavior  of  mild  steel 
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FIGURE  35.  Measured  and  Predicted  Crater  Shapes  for  the  Impact 
of  a 6.35-mm  Sphere  at  60  Degrees  Obliquity  Against  a 6061  Alu- 
minum Alloy  Plate  of  9.53-mm  Thickness  and  at  Various  Impact 
Speeds . 
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Impact  speed,  1.112  km/s 
FIGURE  35.  (Contd.) 
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FIGURE  36.  Measured  and  Predicted  Crater  Shapes  for  Impact  of 
a 6.35-mm  Sphere  at  60  Degrees  Obliquity  Against  a 2024  Alumi- 
num Alloy  Plate  of  6.35-mm  Thickness  and  at  Various  Impact 
Speeds. 
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FIGURE  36.  (Contd.) 
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FIGURE  37.  Measured  and  Predicted  Crater  Shapes  for  Impact  of 
a 6.35-mm  Sphere  at  60  Degrees  Obliquity  Against  a 2024  Alum- 
inum Alloy  Plate  of  3.18-mm  Thickness  and  at  Various  Impact 
Speeds . 


55 


NWC  TP  5844 


Impact  speed,  0.538  km/s 


.Impact  speed,  0.597  km/s 


Impact  speed,  0.644  km/s 
FIGURE  37.  (Contd.) 
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Figures  38  through  44  show  the  experimental  data  obtained  for  firings 
of  6.35-mm  steel  spheres  against  1.47-nim-thick  mild  steel  plates.  In 
Figures  43  and  44  the  solid  line  represents  the  predictions  for  60-degree 
obliqu i ty  using  the  same  value  of  P*  = 4.6  ram  that  was  used  in  Figure  38 
to  fit  the  data  at  normal  incidence.  The  data  in  both  the  ‘ricochet  and 
perforation  ranges  are  widely  discrepant  with  the  prediction  based  on 
this  value  of  P* ^ and  are  seen  to  be  fit  better  by  the  predictions  using 
P*  = 3.0  and  3.8  mm,  respectively,  given  by  the  dashed  lines.  Figures 
41  and  42  show  the  experimental  data  for  impact  obliquities  of  45  degrees, 
and  predictions  with  P*  = 4.6  mm  (solid  lines),  P = 3.5  mm  (dashed  lines, 
ricochet),  and  P*  = 4.4  ram  (dashed  lines,  perforation).  This  shows  that 
the  simple  model  of  ultimate  failure  of  the  plate  using  a constant  value 
of  P*  that  did  well  for  aluminum  is  poor  for  mild  steel. 

There  are  significant  differences  in  the  deformation  and  failure  pro- 
cesses of  mild  steel  and  aluminum  plates  when  impacted  by  steel  spheres. 
These  differences  are  potentially  capable  of  explaining  the  degrees  of 
conformity  to  the  models  proposed  here.  Figures  45  and  46  show  cross 
sections  of  mild  steel  plates  impacted  by  spheres.  In  the  former,  the 
plate  has  been  deformed  but  not  separated  during  a ricochet;  in  the 
latter,  the  plate  has  been  deformed  and  a separation  has  been  formed  dur- 
ing an  embedding  impact.  Figures  47  and  48  show  external  features  of  the 
corresponding  impacts  in  aluminum  targets.  The  mode  of  photography  is 
not  the  same  as  for  the  steel  targets  because  the  overall  failure  process 
is  quite  different.  This  process  consists  of  failures  that  radiate  from 
the  center  of  the  crater  (petaling) . Note  that  the  mild  steel  plate  is 
extensively  stretched  and  thinned  in  the  ricochet  impact,  but  the  alum- 
inum plate  undergoes  the  complete  failure  of  the  petaling  process.  Thus 
it  is  clear  that  the  failure  processes  in  the  mild  steel  target  are  much 
more  prolonged,  and  apparently  are  not  the  simple  deformations  that  are 
modeled  by  the  formula  for  penetration  resistance  nor  the  complete  fail- 
ures modeled  by  the  abrupt  termination  of  resistance. 


CONCLUSIONS 


The  agreement  between  experimental  data  and  predictions  at  oblique 
incidence  shows  that  the  model  of  impact  dynamics  is  reasonably  accurate. 
Comparisons  of  crater  measurements  and  analytical  predictions  show  that 
the  model  of  force  during  ricocheting  impacts  is  a reasonable  approxima- 
tion, but  that  the  approximation  makes  a systematic  error  that  results 
in  overprediction  of  crater  length  at  low  speeds  and  underprediction  at 
higher  speeds.  (There  are  theoretical  bases  for  expecting  a high  tran- 
sient stress  that  is  not  included  in  the  present  model  and  that  could 
account  for  the  discrepancy.)  Comparisons  of  experimental  data  with 
analytical  predictions  of  the  correlation  of  residual  speed  and  obliquity 
with  impact  speed  and  obliquity  show  that  failure  of  the  plate  has  a dis- 
tinct and  identifiable  effect  on  the  correlation  and  that  the  form  of 
modeling  may  be  different  for  different  materials. 
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FIGURE  38.  Data  and  Predictions  of  Exit  Speeds  for  6.35-min 
Spheres  Fired  Against  Mild  Steel  Targets  of  1.47  nun  Thick- 
ness and  at  Zero  Degree  Obliquity. 
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FIGURE  39,  Data  and  Predictions  of  Exit  Speeds 
for  6.  35-tnm  Spheres  Fired  Against  Mild  Steel 
Targets  of  1.A7  mm  Thickness  and  at  30  Degrees 
Obliquity . 


FIGURE  AO.  Data  and  Predictions  of  Exit  Obli- 
quities for  6.35-mm  Spheres  Fired  Against  Mild 
Steel  Targets  of  1.A7  mm  Thickness  at  A5  Degrees 
Obliquity. 
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FIGURE  41.  Data  and  Predictions  of  Exit  Speeds  for  6 . 35- 
nun  Spheres  Fired  Against  Mild  Steel  Targets  of  1.47  mm 
Thickness  and  at  45  Degrees  Obliquity. 


FIGURE  42.  Data  and  Predictions  of  Exit  Obliquities  for 
6. 35-mm  Spheres  Fired  Against  Mild  Steel  Targets  of  1.47 
mm  Thickness  at  45  Degrees  Obliquity. 
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FIGURE  43.  Data  and  Predictions  of  Exit  Speeds 
for  6.35~mm  Spheres  Fired  Against  Mild  Steel 
Targets  of  1.47  mm  Thickness  and  at  60  Degrees 
Obliquity . 


FIGURE  44 . Data  and  Predictions  of  Exit  Obli- 
quities for  6.35-mm  Spheres  Fired  Against  Mild 
Steel  Targets  of  1.47  mm  Thickness  and  at  60 
Degrees  Obliquity. 


61 


NWC  TP  5844 


FIGURE  45.  Cross  Section  of  a Mild  Steel 
Plate  Impacted  by  a Ricocheting  Sphere. 
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FIGURE  46.  Cross  Section  of  a Mild  Steel 
Plate  Impacted  by  an  Embedding  Sphere. 
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FIGURE  47.  Exit  Side  View  of  the  Deforma- 
tion and  Failure  in  a 2024  Aluminum  Alloy 
Plate  Caused  by  a Ricocheting  Sphere. 
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FIGURE  A8.  Exit  Side  View  of  the  Deforma- 
tion and  Failure  in  a 2024  Aluminum  Alloy 
Plate  Caused  by  an  Embedding  Sphere. 
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Specific  features  of  these  correlations,  the  spike  in  the  obliquity 
curves,  and  the  precipitous  drop  in  the  speed  curves  are  associated  with 
failure  processes  in  the  plate.  These  features  are  supported  by  the  ex- 
perimental data,  but  are  not  clearly  obvious  without  the  analytical  pre- 
dictions. 

The  modeling  of  the  force  on  the  sphere  without  failure,  and  the  mod- 
eling of  failure  and  the  effect  of  failure  on  impact  dynamics  are  distinct 
processes  in  this  investigation.  The  recognition  of  these  distinct  parts 
of  impact  dynamics  and  the  role  of  material  behavior  in  each  should  prove 
useful  in  interpreting  further  studies  of  the  correlation  of  impact  pa- 
rameters, configuration  of  the  impacting  system,  and  material  parameters 
with  the  final  outcome  of  impact. 
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NOMENCLATURE 

A Presented  area  of  the  sphere  in  contact  with  the  plate 
B Constant  in  least  squares  fit  of  residual-speed  to  impact-speed 
data 

F Force  on  the  sphere 

Contributions  to  the  force  calculation 
Fp  Component  of  force  that  decelerates  the  sphere 
F]^  Component  of  force  that  deflects  the  sphere 
M Constant  in  the  least  square  fit  of  residual-speed  to  impact- 
speed  data 

P Penetration  of  the  sphere  beyond  the  front  surface  of  the  plate 
Critical  penetration  beyond  which  the  plate  fails 
Radius  of  embedment 
Rq  Radius  of  the  sphere 

5 Distance  along  the  trajectory  of  the  sphere 
T Plate  thickness 

V Speed  of  the  sphere 
Vq  Initial  speed  of  the  sphere 

X Coordinate  of  the  trajectory 

y Coordinate  of  the  trajectory 

a Empirical  constant  in  the  calculation  of  force 

b Empirical  constant  in  the  calculation  of  force 

d Dianeter  of  the  projectile 

h Distance  parameter  in  the  calculation  of  the  relief  of  force  due 
to  plugging 
m Mass  of  the  sphere 

r Radius  of  curvature  of  the  trajectory  of  the  sphere 
t Time 

X Coordinate  in  the  description  of  the  shape  of  the  crater 

y Coordinate  in  the  description  of  the  shape  of  the  crater 

Ps  Angles  used  in  calculation  of  the  presented  area  for  oblique 
impact 

Q Angle  used  to  determine  contact  of  the  plate  with  the  sphere 

a*  Angular  parameter  used  to  locate  failure  of  the  plate 

6 Angular  parameter  used  to  determine  the  surface  of  the  sphere 

that  has  no  resisting  force 

6 Angle  of  obliquity 
Gg  Initial  obliquity  of  the  trajectory 

Angle  used  in  determining  presented  area  after  plate  failure 
J Angle  between  the  force  on  the  sphere  and  the  normal  to  the  plate 

^ Angle  between  the  sphere’s  line  of  flight  and  the  normal  to  the 

contact  surface 
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